Dark-State Polaritons in Single- and Double-A Media 
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We derive the properties of polaritons in single-A and double-A media using a microscopic 
equation-of-motion technique. In each case, the polaritonic dispersion relation and composition 
arise from a matrix eigenvalue problem for arbitrary control field strengths. We show that the 
double-A medium can be used to up- or down-convert single photons while preserving quantum co- 
herence. The existence of a dark-state polariton protects this single-photon four-wave mixing effect 
against incoherent decay of the excited atomic states. The efficiency of this conversion is limited 
mainly by the sample size and the lifetime of the metastable state. 

PACS numbers: 42.50.Gy 



I. INTRODUCTION 

Several years ago, Fleischhauer and Lukin [l[ predicted 
the existence of a stable polaritonic excitation in A-type 
media [Fig. [TJa)] exhibiting electromagnetically-induced 
transparency (EIT) [2]. This excitation, which involves 
a vanishing population of excited states and was there- 
fore dubbed the "dark-state polariton" (DSP), is a coher- 
ent quantum excitation whose evolution is governed by 
a classical control field. This provides a method for ma- 
nipulating single-photon motion, including stopping light 
[3, [4j. The Fleischhauer-Lukin result was derived as a 
perturbation expansion of the field operator equations of 
motion in the strong control field limit. In a subsequent 
work, Juzcliunas and Carmichacl used a Bogoliubov-type 
transformation to diagonalize the model Hamiltonian ex- 
actly, and showed that the DSP can be understood as a 
part of a branch of slow polaritons occurring in systems 
containing a pair of atomic resonances Q. These au- 
thors also emphasized the fact that the photonic part of 
the polariton mixes with atomic excitations possessing 
wavevectors differing by the wavevector of the control 
field. Thus, for instance, it is possible to reverse the 
direction of a polariton wavepacket by switching the di- 
rection of the control field. 

In this paper, we derive the properties of the DSP using 
the Sawada-Brout technique [6|. This can be thought of 
as a simplified version of the method used by Juzeliunas 
and Carmichael, and we shall see how its results reduce 
to those of Fleischhauer and Lukin near resonance, which 
was not demonstrated in Ref. [5j. We then extend the 
analysis to a double-A medium [Fig. [T^b)], which con- 
tains a DSP consisting of low-lying atomic excitations 
and photon states of two different frequencies @, [H ■ In 
both single- and double-A systems, the DSP is protected 
against incoherent decay processes acting on the excited 
states, because it contains a vanishing population of these 
states. It has previously been shown that double-A me- 
dia can efficiently upconvert classical probe beams [gLfToj] . 
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FIG. 1: (a) 3-level A-type medium, (b) Double-A medium. 



and a related four-wave mixing scheme has already been 
used in such systems to generate correlated photon pairs 

ELI GJ, H G3- Here > we p° int ° ut that the DSp could 

be exploited to perform single-photon frequency conver- 
sion in a manner that preserves quantum coherence. Un- 
like semiclassical analyses in which the electromagnetic 
field is treated classically @, [l(| , this theory applies to 
the single-photon regime. It may thus have applications 
in quantum information processing, such as for down- 
converting a member of an entangled photon pair to a 
frequency suitable for transmission over a telecommuni- 
cations fiber. Unlike parametric conversion schemes ex- 
ploiting optical nonlincarities, the relevant photons are 
up- or down-converted individually, instead of being split 
or recombined; the additional momentum and energy are 
supplied by the control fields. 



II. SINGLE-A SYSTEM 

We begin by considering an Y-atom gas with a single-A 
level structure, shown in Fig.QJa). The ground, excited, 
and metastable atomic states are respectively denoted by 
I b), I a), and |c), and their corresponding energies by fkOb, 
huj a , and fwj c . The atomic Hamiltonian is 



H = hJ2(u a ar 



hi 



(1) 



with the sum performed over all atomic positions r. Here, 

WrR (2) 
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denotes a transition operator for the atom at position 
r. We also define Fourier-transformed operators, ai b = 



2 



N 1 / 2 J2 r &r b e%kr etc - The photon Hamiltonian is 

Hi = ?ic|fc|a^a fe , (3) 
fc 

where aj. and a k are photon creation and destruction 
operators. The photons interact with the ab transition 
through the minimal-coupling Hamiltonian 
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h.c. 



(4) 



The coupling constant is g ~ Vy / 2nNuj a b/hV, where V 
is the dipole moment of the ab transition, u) a f, = Lu a — u) b , 
and V is the cavity volume. For notational simplicity, 
we have used the rotating-wave approximation. Finally, 
we include a classical control field with strength (Rabi 
frequency) f2, frequency ujl ~ ^ac, and wavevector k^'- 



H 3 (t) 



-hfl e 



^ e ikLr < c - 



h.c. 



(5) 



Here, we have again discarded counter-rotating terms. 
We neglect the coupling between the photons and the ac 
transition, which is negligible compared to the effects of 
the control field, and the coupling between the control 
field and the ab transition, which is off resonance. 

The time dependence in ([5]) can be removed by defining 

H L = U L {t)H{t)ul(t) + Tiuj l J2 < c , (6) 

r 

where H(t) = H H h H 3 (t), and 

U L (t) = exp ~iuj L tJ2< C ■ (?) 



The Schrodinger equation Hit) \ip(t)) = ifidt 
then rewritten as 



in- [u L (t)\m)\ 



h l [u L {t) m))} 



can be 



(8) 



Thus, we can extract solutions to the Schrodinger equa- 
tion from the energy eigenstates of the time-independent 
Hamiltonian Hl. To obtain these, we look for a polariton 
excitation operator A^ such that 



[H L ,A^] = TiluA ] 



(9) 



The polariton is long-lived provided the omitted terms 
are negligible [f| . If the initial state of the system is its 
(zero photon) ground state, should be a mixture of a T , 
a ab , and a cb . Below, we list the commutation relations of 
these three operators with Hl ■ We have removed terms 
involving a ba , a aa , a ca , and a^ since these operators 
give zero when acting on the ground state, this intro- 
duces no additional error for single-polariton excitations. 
Similarly, we have replaced a k b with \/N5ko- Thus, 



[Hl,o?] 
[HlA] 



h(uj cb + uj L )cr c k h -k L 



hc\k\a\ — hga 
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- hg*4 ( 10 ) 
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FIG. 2: Polaritonic dispersion curve for ujl = Woe- The solid 
lines show the exact polariton solutions given by Eq. (|14p ; the 
dashed line showss the Fleischhauer-Lukin solution, Eq. (|15l) . 



Let us now look for excitation operators of the form 



.4! 



il ab 



l „cb 
ik "k-kr 



+ Kk<, (13) 

where the band index n enumerates the different polari- 
ton species. The c-numbers <p 3 nk are determined by in- 
serting ||T5D into © and using QI5|)-|(T2J). This gives three 
self-consistency equations that can be written as 
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(14) 
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The form of the effective Hamiltonian in fTJJ is familiar 
from semiclassical analyses of EIT. Fig. [5] shows the band- 
structure in the absence of loss, similar to the one given 
in Ref. ||. For simplicity, let us assume that W£ = oj ac . 
The asymptotic eigenfrequencies far from resonance are 
c\k\ and the eigenvalues of the upper-left 2x2 submatrix 
in the effective Hamiltonian, in this case uj ab ± f2. Ex- 
actly at resonance (|fe| = uj ab /c), there is an eigenvector 
oc [0, 1, — Cl/g], and for slightly detuned k this eigenvector 
continues into ones where the a ah component is nonzero 
but small. These solutions — "dark-state polaritons" — 
are thus insensitive to incoherent decay processes acting 
on | a). The stability of the exactly-resonant DSP is lim- 
ited only by the lifetime of the metastable state |c) , which 
we shall assume to be longer than the time-scale of any 
relevant experiment. For off-resonant DSPs, the decay 
rate is only quadratic in the detuning: upon replacing 
Lu a with uj a — iT a in (|14[) . one finds that the imaginary 
part acquired by uj k is ~ r a |A/f2| 2 (for 3> g), where 
A = c\k\ — u> ab . The other two polariton branches corre- 
spond to "bright" polaritons that contain significant \a) 
population and are thus strongly affected by losses. As 
in Ref. [l|, we neglect Langevin noise effects, which do 
not influence the adiabatic evolution of the DSPs. 

Expanding around u = uo ab yields a limiting solution 
for the DSPs: 
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(c\k\ - UJab) 
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(15) 
(16) 
(17) 
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Eq. (|15H17[) agree with the solution derived by Fleis- 
chhauer and Lukin using a perturbation expansion in 
I/O In our formalism, the fact that decreasing O 
reduces the polaritonic group velocity can be intuitively 
understood as the result of "squeezing" the bandwidth 
of the middle polariton band. An interesting property of 
the DSP solution is that it does not depend on the ener- 
gies of the underlying A system, only the detuning of the 
control field and the coupling parameters g and Q. 

Finally, we can extract the solutions to the original 
Schrodinger equation using (J5J). For a polariton with 
quantum numbers (n, fc), the state at time t is 



-0Vf +<& k j" L 'of-k L +€k4 |0). (18) 



The a cb component in (fl8| has a different frequency and 
wavevector compared to the rest of the polariton. This 
property does not, however, destabilize the polariton: in 
a wavepacket constructed of a superposition of DSPs, 
the photonic and a cb components possess different phase 
factors but share a single envelope. 

The preceding derivation holds regardless of the angle 
between the input photon and the control beam. The 
direction of kh only enters into the choice of excitation 
operator <J k b _ kL occurring in the polariton operator (p~3|) . 
and plays no role in the eigenproblem (fT4|) that yields the 
state amplitudes and polariton energy. 

By switching between two non-collinear control beams, 
it is possible to coherently rotate the photon wavevector, 
by an angle of up to 2sin~ 1 ((x> ac /ti; a k), where the plane 
of rotation is specified by the polarization of the con- 
trol field. A special case of this has been discussed by 
Juzeliiinas and Carmichacl: when cut ss uj c , one can co- 
herently backscatter the photon by inserting a photon 



with k ~ kjj, which mixes with a a excitation with 
wavevector k — k^ ~ 0, and switching the control field 
to —kjj. The a bc excitation then mixes into a photon of 
wavevector k — 2k^ ~ — k 



III. DOUBLE-A SYSTEM 

Suppose we add a second excited state, \d), as shown 
in Fig. Hfb). A second control beam couples \d) to \c), 
and for simplicity we assume that the two control beams 
have parallel polarization vectors. The d <-*■ a transi- 
tion is assumed to be forbidden. One of the reasons this 
"double-A" system is interesting is that it can be used 
to upconvert or downconvert probe beams, as experi- 
mentally demonstrated by Merriam et. al. [10( and other 
groups. It can be shown, using the Fleischhauer-Lukin 
formalism, that this type of level structure supports a 
DSP [![. As we shall see, this DSP arises naturally from 
the present method as a 5 x 5 generalizion of (TTi|) . 

The Hamiltonian, H'(t), contains four new terms. The 
first, Y^r w d a r d i gives the energy of the \d) states. Next, 
we introduce a second photon field with operators bt and 
bk, and Hamiltonian Tic\k\b* k bk- (There is really only 
one photon field, but this trick is permissible since the 
atom-photon coupling becomes negligible far away from 
the EIT resonances.) Finally, we add interaction terms 
analogous to Q and ([5]), with \d), bk, g', and tt' replacing 
| a), a*,, g, and fl respectively. 

The control field interaction Hamiltonian (J5]) and its 
analog for the dc transition oscillate at different frequen- 
cies, so the transformation ©-([7]), which works by ro- 
tating |c), cannot eliminate the time dependence. We 
can overcome this difficulty with a transformation that 
instead rotates the \a), \d), and photonic states. Let 



h' l = u' L {t)H'{t)u'j} - hcj L (J2 k a Uk + 

where H'(t) is our new Hamiltonian, and 

u' L = cx P [iu> L t (J2 k a Uk + J2< a ) +^ (E fe 5 !^ + E r a ' 



d<i 



(19) 



(20) 



This once again allows us to write the Schrodinger equation as ihdt [U' L (t) \ip(t))] = H' L [U' L (t) \tp(t))], where H' L is 
time-independent. We look for excitation operators for H' L of the form 



Al 



a b 



l2 „db 
J nk a k+k' 



b nk a k 



4 t 
nk a k+k L 



(21) 



The self-consistency equations for the parameters 

Hi = h 



b J nk take the same matrix form as (|14[) . with effective Hamiltonian 
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Udb — u' L 
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-k L \- LU L 
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c\k- 


Yk' L 



(22) 



The polariton created by (|2"Tj) is a valid excitation be- 



cause, as in the single-A case, no extra non-negligible 
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Wcb/C 



This equation is exactly satisfied by (|2ip . for which 
lu = c(k + k^) and lu' — c(k + k' L ). The physical mean- 
ing of (f2"8")) is particularly easy to deduce in the present 
theory: in the single-photon limit, the stationary state 
corresponds to a polaritonic solution of the form (|2ip . 
for which the photonic components cannot take on arbi- 
trary frequencies because they are coherently mixed via 
the atomic excitation a^ b . The amplitude-matching con- 
dition for the semiclassical dark state is 



FIG. 3: Polaritonic dispersion curve for the double- A medium. 
The solid lines show the exact polariton solutions given by 
Eq. (I22|l ; the dashed line shows the linearized solution given 
by Eq. (|23[) . The horizontal asymptotes occur at ui c b and 

u cb ± VN 2 + |n'| 2 . 



terms are generated by commutating this operator with 
the Hamiltonian. Let us now assume that the control 



U)„ 



and uj'r 



For 



fields are resonant, i.e. u>l 
\k+k^\ = Uab/c and |fc+fc^| = ujdb/c, the effective Hamil- 
tonian (|22p has an eigenvector tx [0, 0, 1, — Q/g, — ft'/g'}. 
The first two components of this eigenvector, correspond- 
ing to the two excited states, are identically zero, so this 
represents a DSP consisting of crjjr excitations and pho- 
tons with wavevectors k + k^ and k + k' L . It can be shown 
that no other linearly independent eigenvector with this 
property exists, so there is only one such DSP solution. 
The linearized DSP solution, analogous to (fT5")) - (fl"T)) . is 

\n/g\*6k+\n>/gT6k> 

(24) 

(25) 

(26) 
(27) 

where 8k = \k + k^\ — uj a b/c and 5k 1 = \k + k' L \ — LUdb/c. 
The polaritonic bandstructure, in the absence of loss, is 
shown in Fig. [3l 

These results can be shown to be consistent with 
the single-photon limit of a semiclassical analysis of the 
double- A medium given by Korsunsky and Kosachiov Q . 
In this one-dimensional model, where the electromag- 
netic field is treated classically, the Heisenberg equations 
of motion for the atomic system possesses a stationary 
"dark state" solution that is decoupled from the electro- 
magnetic field and is stable against spontaneous emis- 
sion. This dark state exists only if the background field 
(consisting of two probe beams and two control beams) 
obeys certain frequency, amplitude, and phase matching 
conditions. The frequency-matching condition is 



^nk 


= Ucb + — 


f \n/ g \i + \n'/g'\2 
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\n'/g'\ 2 {8k-Sk') 
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\g\ 2 i + 


\Sl/g\* + |fi'/<?1 2 
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n' 8k' 
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= -(fi'/s') 







lu ■ 



^cb, 



(28) 



where lu and lu' are the respective frequencies of the 
probe beams resonant with the ab and db transitions. 



VE 



VE' 

n' 



(29) 



where V and V' are the dipole moments for the ab and 
db transitions, and E and E' are the electric field ampli- 
tudes of the associated probe beams. The electric field 
amplitudes can be related to the quantum mechanical 
photon amplitudes <^> 4 and <fi 5 by 



E ^ y / 2TThuJab/V cj) 4 

E' «-» ^/2Trhuj db /V<f> 5 , 



(30) 
(31) 



which can be verified by computing the expectation value 
(|_E| 2 ) produced by each photon creation operator. With 
this identification, the linearized DSP amplitudes (j2"6")l 
and l|27p satisfy (|2^|) . The third condition derived by 
Korsunsky and Kosachiov, which relates the phases of 
the four beams, is also satisfied by the DSP because, as 
shown by (f2"6")) and (|2"T|) . the phases of the probe beams 
are locked to those of the control beams and fi'. 

The dark state studied by Korsunsky and Kosachiov is 
a pure state of the atomic system, reflecting the fact that 
the electromagnetic field is treated classically @ . In con- 
trast, the present model takes into account the coherent 
mixing between the quantum state of the probe field and 
the quantum state of the atomic medium: performing a 
partial trace of the DSP over the photonic Hilbert space 
yields a mixed atomic state. This mixing becomes impor- 
tant at the single-photon level, which is also potentially 
the regime of interest for quantum information process- 
ing. In the following section, we will examine how this 
mixing can be used to convert between the two photonic 
components of the double-A DSP. 



IV. FREQUENCY CONVERSION 

For a single-A medium with a resonant control beam, 
inserting a photon with wavevector ko, resonant with the 
ab transition, gives rise to a DSP whose group velocity 
points in the same direction, independent of the direction 
of the control beam. This freedom to choose the direc- 
tion of the control beam disappears in the double-A case. 
Here, an incident photon fco mixes with another photon 
with wavevector k\ = ko — k^ + k' L . Assuming both con- 
trol beams are tuned to resonance, the resulting state 
overlaps with a DSP only if \ki\ ~ LOdb/c. Furthermore, 



5 



the group velocity of the DSP is, from (|23[) . 

\n/g\* fc + \W/g'\ 2 h 



1 + \n/g\* + \W/gf 



(32) 



where fco = fco/|^o| and k± = ki/\ki\. Therefore, a choice 
of fco and fci determines the directions of the two con- 
trol beams (or, more generally, choosing any two of these 
directions determines the other two). As an aside, we 
note that the beam matching conditions forbid the choice 
fco = —ki, which would imply the possibility of a station- 
ary wavepacket with nonzero control beams; however, if 
fco and ki are nearly antiparallel, (I32|) predicts that the 
control beam strengths can be tuned to produce a low 
group velocity. 

In order to illustrate the mixing between the two pho- 
tonic components in the DSP, let us fall back on the 
"trivial" one-dimensional case where all wavevectors are 
parallel, which satisfies the above beam matching condi- 
tions. Suppose we inject the photon fc at t = 0, so that 
the quantum state is 



|^(0)) 



10) 



(33) 



where k = fco — fcz,. Without losses, the state at time t is 



At (e- WWR ) 54 &Uj°> + 



(34) 



where the matrix TC is defined in (f2"2"|) and the omitted 
terms are the other polariton components. The result, 
shown in Fig. [H is an oscillating upconversion amplitude 
| (0| bk+k' L I that can approach 100%. The effects 

of incoherent excited state decay, which can be modeled 
by replacing uj a with ui a — iT a and u>d with oj a — i^d in 
(|22)l . are also shown in Fig. |4j Although the DSPs are 
protected against decay, damping still occurs because the 
incident photon generates a non- vanishing population of 
bright polaritons. When these exit the system (typically 
as off-axis photons), only the DSP remains. 

A more efficient example of single-photon frequency 
conversion can be obtained by going from momentum 
space to real space and studying the behavior of polariton 
wavepackets. Let us define c-number fields $j — &j(r, t) 
such that 



+$ 3 of + $4 e- ikhT a\ + $ 5 e-**WtJ . (35) 

Inserting (|35[) into the Schrodinger equation and using 
(|22p . we obtain a Schrodinger wave equation 



9$ 



(36) 



If k is chosen such that |« + &l| = ^ab/c and \k + k' L \ = 
uidb/c, then the DSP corresponds to values of $>j that 




t(10" 3 s) 

FIG. 4: Numerical solutions of | (0| b k _ kL+k i \ip(t)) | against 
t, where \tl>(t)) is the quantum state at time t after inserting 
a photon at with k = CJ a b/c- Here, cu c b = 10 4 cm -1 , \g\ = 
\g'\ = 0.1 cm -1 , and |fi| = \Q'\ = 1cm -1 . (a) No excited 
state decay, r a = T d = 0. (b) r a = F d = 0.02 cm" 1 . 





-5 0.6 
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FIG. 5: Photon frequency conversion. The a-photon ampli- 
tude \4>4,(z, t)\ (dashed line) and &-photon amplitude |0s(.z,i)| 
(solid line) are plotted at three instants. The abscissa is 
z/ cosd, where kh • z = cos#. The second control beam is 
pointed that k' L ■ z — 0.9 cos 9. Both control beams are c.w. 
The effective thickness of the double-A medium, zo = 60 cm, 
is indicated. Within the medium, \g\ = \g'\ = 0.1 cm -1 , F a = 
T d = 0.2g, |OJ = 1cm" 1 , and |fi'| = 3{1 + tanh[4(z/ cos 6 - 
0.5)/zo]}cm _1 . Outside the sample lies vacuum. The ampli- 
tudes are computed by integrating Eq. (|36[) numerically for 
the initial conditions in Eq. (|38|l . where f3 = 0.8 m -2 and 
Z2 — —50 cm. 



are constant in space. For a wavepacket centered around 
k with bandwidth <§; w a fc, uidb (he., spatial width much 
longer than the optical wavelength, which is the usual 
slowly- varying envelope approximation), TL takes on the 
intuitive local form 



H(r,t) « h 
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n* 
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n g 

n' o 

o o 

— icfco 







9' 




— zcfci 



(37) 



As in the single-A case, the evolution of the polaritonic 
envelope is independent of the underlying double-A fre- 
quencies, except through the coupling parameters g, g' , 
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FIG. 6: Frequency conversion efficiency, parameterized by 
the converted photon amplitude $ 5 normalized to the in- 
put photon amplitude, as a function of the sample length 
zq. Within the sample, the control beam Q! varies as = 
3{1 + tanh[4(z/ cos 6 — 0.5) /zq]} cm -1 . Curves for incoherent 
decay rates r a = Fd = 0.2g, OAg, and 0.6g are shown. All 
other parameters are the same as in Fig. [5] 



fl, and n'. We again emphasize that this result is not 
perturbative; it holds for arbitrary values of f2 and fl', 
and depends only on the fact that the wavepacket is suf- 
ficiently broad. Generally, the coupling parameters can 
vary (smoothly) in space; for instance, a variation in fi 
or f2' could be accomplished using a c.w. control beam 
with a non-uniform cross-sectional intensity profile. Such 
variations can be used to "adiabatically" transfer one 
photon population to another within a propagating DSP 
wavepacket, substantially improving the efficiency of the 
conversion process compared to the previous example. 

Let us consider an effectively one-dimensional experi- 
mental setup where all relevant spatial variations occur 
in the z direction. In particular, we must assume that the 
x and y edges are far enough away that boundary effects 
(which appear when the beams are not all collinear) are 
negligible. The incident envelope field &j(z, t = 0) is 



<E>4 = exp 



COS ( 



Z-2 



(38) 



$1 = $ 2 = $3 = *5 = 0. 



Outside the sample (z < or z > zo), all coupling pa- 
rameters are zero. Within the sample (0 < z < zo), the 
functional forms of Sl(z) and Cl'(z) are chosen so that 
1^1 ~ \d\ ~ l^'l near the entrance of the sample, which 
ensures that the DSP is dominated by the input photon; 
whereas |0'| > > \g\ near the exit, which ensures that 
the DSP is dominated by the converted photon. The re- 
sult is shown in Fig. [5) For the given parameters, the con- 



verted photon amplitude is ~ 0.9 times the incident am- 
plitude. The efficiency is limited by the available length 
of the double-A medium. As shown in Fig. a longer- 
sample allows the f2' field to be varied more gently, gen- 
crating fewer bright state polaritons and increasing the 
conversion efficiency. 

V. CONCLUSIONS 



In this paper, we have presented an analysis of single- 
and double-A EIT systems based on a microscopic 
equation-of-motion technique. Within this formalism, 
the presence of a DSP corresponds to the existence of 
special eigenvectors of an effective Hamiltonian matrix, 
in which the entries corresponding to rapidly decaying ex- 
citations are identically zero, regardless of the strength 
of the control fields. The ability of the double-A sys- 
tem to efficiently upconvert and downconvert photons, 
previously established in semiclassical four-wave mixing 
studies [1, [l(| , is retained in the coherent single-photon 
limit due to the existence of the DSP. The analysis can be 
further generalized to multi-A systems, where one finds 
additional polaritonic bands similar to those in Fig. [3J 
with exactly one family of DSP solutions possessing van- 
ishing excited state populations. 

Throughout this paper, we have restricted our atten- 
tions to the single-polariton sector of the theory, which 
is valid only if the polaritons are much more dilute than 
the underlying atomic medium. The polariton operators 
(fT"3f and (|33[) do not obey exact bosonic commutation re- 
lations, since the a operators are not bosonic operators; 
however, the corrections to the commutator vanish as 
O(MfN), where M is the number of atoms excited [l~5|. 
This condition is satisfied, for instance, in the experi- 
ments of Merriam et. ai, where M/N ~ 1(T 3 [10]. The 
single-polariton sector has the advantage that the quan- 
tum state of the system can be expressed in terms of a 
simple wave equation, as in (|36|) . Thus, once the a oper- 
ators have been used to derive the effective Hamiltonian, 
the additional structure given by their non-trivial com- 
muation relations disappears from the theory. Should one 
wish to study the limit where M becomes comparable to 
N, this structure will have to be taken into account. 
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